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1 Introduction 



Quantum R matrices, solutions of the Yang-Baxter equation, are interpreted as diffusion amplitudes 
for two-body interactions of particle type eigenstates in integrable two-dimensional field theories. 
The Yang-Baxter equation then ensures consistent factorisability of the three-body amplitudes in 
terms of two-body ones. 

When considering integrable field theories with non trivial boundary effects such as theories on 
a half-line, one needs to introduce a new object describing reflection processes on the boundary. 
Integrability is preserved provided that the two-body exchange R matrix and the one-body reflection 
K matrix obey a quartic consistency condition [H El El H] 

R 12 (u - v) K^u) R 12 (u + v) K 2 (v) = K 2 (v) R l2 (u + v) K^u) R 12 (u - v) . (1.1) 

Let us remark that the reflection equation also appears when considering generalisation of ZF algebras 
(i.e. Zamolodchikov-Faddeev algebras allowing the presence of a boundary. These generalised 
ZF algebras ensure the total scattering matrix of the model to be unitary. In this approach, the com- 
mutation relations of ZF generators implement an operator b(u) which obeys the reflection equation. 
The K matrix would then be a representation of this boundary operator b [H] . 

Another point of view has recently been presented using a universal construction of reflection 
algebras as twists of quantum algebras [ZHHj- The exchange equation (jl.lj) would then describe a 
trivial representation of the quantum generators, although it is not clear to us that this interpretation 
is valid for spectral parameter dependent solutions. Furthermore, solutions with quantum degrees of 
freedom on the boundary have also been obtained in 

These reflection equations, or boundary Yang-Baxter equations, have recently drawn attention 
and systematic ways of computing some solutions for given R matrices were derived e.g. in jTT H lll j lTS*]. 
recovering and extending previous results derived for instance in f3J El El E3- The considered 
cases were related to trigonometric R matrices for all spins 

p], and to p2] and A [ n ] vector 

representations ^2], yielding a wealth of new solutions. 

In this paper, we will consider the case of R matrices corresponding to vector representations 
of Yangians and super- Yangians, i.e. rational R matrices constructed in [TBI HH]- The complete 
classification for gl(n) in vector representation was given in j20]: it appears that generically any 
solution is conjugated (by a constant matrix) to a diagonal solution. [21] then derived a series of 
solutions for so(m) and sp(n) based upon an ansatz proposed by Cherednik pp. We present here, 
for the Lie (super) algebra series so(m), sp(n) and osp(m\n), a classification of purely diagonal, 
purely antidiagonal and mixed diagonal/antidiagonal solutions by directly solving the boundary 
equation (jl.lj) for a one- dimensional boundary quantum space. Once the diagonal case is exhausted 
(corresponding to "flavour-preserving" reflection matrices), the most natural extension to look for 
indeed consists of reflection matrices which preserve pairs of conjugate states (according to (j2.2jl ). 
Particular solutions have already been derived for the so{m) algebra [211 E2j and the spin) algebra 
j21j . They can all be identified with particular diagonal solutions of our classification. 

We must point out that there exists in fact another notion of reflection equation, which is related 
to the definition of twisted Yangians [2S1E11, and arises also in the theory of coideal algebras described 
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in [23 12H1 • This equation reads 



R 12 (u - v) K^u) R%{-u - v) K 2 (v) = K 2 (v) R%{-u - v) K^u) R l2 {u - v) 



(1.2) 



where the transposition t is defined below (see definition 12. 2|) . As pointed out in [THj, this equation is 
actually the same as when R{u) is the R matrix of Yangians of type so(n), sp(n) or osp(m\n). 

One essential purpose in establishing such a classification of reflection matrices is to use them in 
order to construct and eventually to solve spin chain models with a variety of boundary conditions. 
We present here explicit resolutions, using the analytical Bethe Ansatz method, of the so(n) and sp(n) 
open spin chains with boundary conditions determined by the diagonal solutions of the reflection 
equations. In particular, we derive explicit formulae for the eigenvalues of the low-lying excitations of 
hole type. We obtain in each case the full scattering matrix without ambiguities (e.g. CDD factors) 
including bulk and boundary interactions. This result has in fact a very general character, not 
limited to the particular models considered here. This actually provides us with bulk and boundary 
S matrices relevant for integrable field theories with non trivial boundary conditions, for which they 
represent "universal" S matrices. We expect, in analogy with the bulk case, that such boundary S 
matrices should correspond to sp(n) or so(m) Gross-Neveu model with certain boundaries. For this 
particular model, the study of boundary conditions associated with the different K matrices deserves 
further attention but goes beyond the scope of the present paper. 

2 Generalities 

Let gl(m\n) be the Z 2 -graded algebra of (m + n) x (m + n) matrices Xij and 9 = ±1. The Z 2 - 
gradation is defined by (— 1)^ = 9o if 1 < i < m and (— 1)^ = — #0 if Tn + 1 < i < m + n. In the 
following, we will always assume that n is even. 

Definition 2.1 For each index i, we introduce a sign 9i 




for 
for 




1 < i < m + 



n 



2 



i < m + n 



(2.1) 



and a conjugate index 1 




for 1 < % < m 

for m + l<i<m + n 



(2.2) 



In particular 9i9 T = 9q( 



Definition 2.2 For A 



. A l i Eij, we define the transposition t by 




(2.3) 



R satisfies (A 1 ) 1 = A and, for C-valued bosonic matrices, (AB) 1 = B t A t . 
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As usual Eij denotes the elementary matrix with entry 1 in row i and column j and zero elsewhere. 
We shall use a graded tensor product, i.e. such that, for a, b, c and d with definite gradings, 
(a <g> 6)(c® d) = (-l)MWac ® bd. 

Definition 2.3 Let P be the (super)permutation operator (i.e. X 2 % = PX\ 2 P) 

m+n 

P=Y d (-l)®E ii ®E ji (2.4) 

*>j=i 



and Zei 



We define the R matrix 



Q=J2 (-l) m 0i8jE Tl ® £7^ = P* 1 . (2.5) 
-R(m) = I + — — — (2.6) 

U U + K 



with Ik = (m — n — 2)8q. 

The R matrix (|2.6j) satisfies the super Yang-Baxter equation 

Ru(u) Ri 3 (u + v) R 23 (v) = R 23 (v) R 13 (u + v) R 12 {u) (2.7) 

where the graded tensor product is understood. 
The operators P and Q satisfy 

P 2 = 1, PQ = QP = 6 Q, and Q 2 = 6 (m - n)Q . (2.8) 

The R matrix ()2.6j) is known to yield the osp(m\n) Yangian ^H], arid leads to the non-super or- 
thogonal (taking n = 0, 6 = 1) and symplectic (m = 0, 6q = — 1) Yangians. For obvious rea- 
sons, we will call labels as orthogonal (resp. symplectic), indices i which satisfy 1 < i < m (resp. 
m + l<i<m + n). 

Although we will restrict ourselves to the diagonal, antidiagonal and mixed cases for the reflection 
matrices K, the following lemma can be used to get more general solutions. 

Lemma 2.4 Let K{u) be a solution of the reflection equation M.l)) and U such that UU l = 1 be a 
(constant) matrix of the orthogonal, symplectic or ortho symplectic (super)group (depending upon the 
choice of R). Then K l (u), UK{u) V 1 and U^iu) U 1 are also solutions of M.l\) . 



The proof is straighforward, using the invariance of the R matrix under conjugation by U and 

E>tlt2 D 

"12 — -"-12- 



3 



3 Diagonal solutions of the reflection equation 

In this section we consider invertible diagonal solutions for K(u), i.e. 

K(u) = diag (fci(u), • • • , k m (u) ; k m+1 (u), ■■■ , k m+n (u)^j (3.1) 

where the semicolon emphasises the splitting between orthogonal and symplectic indices. Here 
the ki(u) are supposed to be analytic C-functions of u, i.e. the boundary quantum space is one- 
dimensional. 

Proposition 3.1 There are three families of generic diagonal solutions and two particular cases 

Dl: Solutions of sl(m + n) type, with one free parameter, form even 

1 cu ui n 

fci(u) = l, Hu) = - , Vi G {1, to + 1, to + — } (3.2) 

1 — cm 2 2 

27ns solution has no extension to odd to. 

This solution is obviously invariant under the action of 5X(y||). 
D2: Solutions with three different values of ki{u), depending on one free parameter 

u( u) - 1 + c ^ u ( u )- 1 + CmU 

1 V / j rvin\a) — 

1 — C\U 1 — c m u 

kj( u ) = 1 7^ 1, m 

where (k — )cic m + C\ + c m = (3.3) 
27ms solution does not hold for m = 0,1. 

27ie moduli space of this solution is invariant under the action of OSP(m — 2|n) x 5*0(2). 
.Dc?.- Solutions without any free continuous parameter (but with two integer parameters) 
ki(u) = ki{u) = 1 Vi G {1, mi; to + 1, to + rii} 

fcj(u) = kAu) = V? G {toi + 1, to — toi; to + m + 1, to + n — ml 

1 — cu 

2 

where c = — (3.4) 

k - 9 (2m 1 - 2m - 1) 

27ws solution is invariant under the action ofOSP(2mi\2ni) x OSP{m — 2mi\n — 2n\). 
DJf.: In the particular case of so(4) 7 the solution takes the more general form: 

ts< \ A- fi l + c 2U 1+C 3 M l + C 2 ul + C 3 u\ 

K[u) = diag 1 , , , 3.5 

\ 1 — C 2 U 1 — C3M 1 — C 2 U 1 — C 3 M / 

This solution contains the three generic solutions Dl (c 2 c 3 — 0), D2 (c 2 + C3 = 0) and D3 
(c 2 = c 3 = 00). 



D5: In the particular case of so(2), any function-valued diagonal matrix is solution. 



In each case D1-D4, infinite value of the parameter c is allowed. It corresponds to the constant 
value i±2 = -1. 

l—cu 

These solutions are given up to a global normalisation and a relabelling of the indices, provided 
it preserves the orthogonal/ symplectic splitting and all the sets of conjugate indices {i,i}. 

Let us remind that the solutions for so(m) and spin) algebras are obtained by setting n = or 
m = respectively. 

Solutions found in [22] for so(m) algebras (m > 2) are identified after a suitable change of basis 
with the sets Dl and D2. Solutions found in [21] for so(m) and sp(n) algebras are identified with 
the set D3 and limiting cases of Dl. The solution D3 with the special value m\ = (m — l)/2 (with 
n = 0) corresponds to the only non trivial rational limit of the set of trigonometric solutions given 
in [2B!. 

Remark that the generic solutions Dl, D2, D3 are associated with symmetric superspaces (as they 
were introduced in [29J, in a different context) based on OSP(m\n), namely OSP(m\n)/ '5L(y|^), 
OSP(m\n)/OSP(m - 2\n) x SO(2), and OSP(m\n)/OSP(m - 2m 1 \n - 2n x ) x OSP(2m 1 \2n 1 ), 
respectively. In the case of Lie algebras, we recover the usual symmetric spaces (based on orthogonal 
and symplectic algebras). 

Proof: The projection of the reflection equation (jl.lj) on Ey <8> Eji for i 7^ j, J reads 

ki{u)ki(v) H kj(u)ki(v) = kj{u)kj(v) H kj(v)ki{u) , (3.6) 



u + v u—v u + v u—v 

the solution of which is given by 

ki(u) _ 1 + Ciju 



kj(u) 1 — CijU 



i^jj- (3.7) 



The cases Qj = 0, 00 correspond to the constant ratios |^ = ±1. 

For convenience, we will introduce Fij(u) = (well-defined since K(u) is supposed to be invertible 
for generic u). When i 7^ j, J, is then given by ()3.7j) for some qj and obviously Cji = — Cjj. When 
defined, since Fij(u)Fjk(u) = Fik(u), the parameters c must also satisfy 

for iy^j,], j 7^ k, k and k 7^ i,i . (3.8) 

CijCjkCki 

To any solution for K(u), we associate a partition of {1, . . . , m + n} where the classes are defined by 
i = j <^ ki — kj (•<=> c^ = 0). The constraints ()3.8j) are sufficient to conclude that the partition 
associated to any solution has at most three different classes, except when m + n = 4, where it can in 
principle have four classes. Note that in the case of sl(m\n) where the constraints ()3.8)1 hold without 
any restriction on the indices, K(u) is always built with at most two different functions. 
More precisely, if m + n 7^ 4, ()3.8|) implies that the partition of indices is constituted either of two 
subsets, or of three subsets. In the last case, it can only be of the form 
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(D2) {i}, {i} and {1,..., m + n} \{i,i} 

where i is an orthogonal index, i.e. iG {1, . . . , m}. 

Projecting the reflection equation (jl.lj) on E^® E^, one gets (for i ^ after taking u — > v 
k {E Ti {u) - Fjjiu)) + 6 {Fjiiu) - Fijiu)) + (2u + k) (F^u) - F- xj {u)) 

= (3.9) 

i 

This equation show in particular that — 1 F^ =1. A recursion allows us to write the 
two-subset partitions as either 

(Dl) I = {z*i, i 2 , • • • } and / = ■ ■ ■}, i.e. i & I I 

or 

(D3) / = {ii,H,i 2 ,h ■■■} and J = {jijij 2 ,j2, ■■■}, i.e. iE I OiE I 

Knowing the possible forms of the partition, one can evaluate the sum in ()3.9j) . This equation involves 
at most two different functions F mn . It provides the constraints on the parameters c mn . A global 
check ensures that all the remaining projections do not lead to new constraints. 
The cases m + n = 4 are solved by direct computation. Only the so(4) case eventually exhibits a 
more general solution (D4). 

Finally, the case of so(2) is special. In particular, R(u) appears to be diagonal. A direct computation 
shows that all function-valued diagonal matrix is solution. ■ 
Note that when k = (as in the case of so(2)), the spectral parameter can be rescaled before 
taking the limit k — *■ 0, and the corresponding R matrix does not involve the identity anymore |3Uj . 
This is the R matrix used in [22] for so(2). 

4 Antidiagonal solutions of the reflection equation 

We now look for invertible solutions of the "antidiagonal" form 

m+n 

K{u) = Y d l i {u)E ti . (4.1) 

i=i 

Proposition 4.1 Solutions of the reflection equation with antidiagonal terms exist only in the pure 
so(2m) or sp(2n) cases. They are constant and only restricted by the set of constraints 

= 1 , Vi. (4.2) 

In the special case of so(2), any function-valued antidiagonal matrix is solution. 

The proof of this proposition will be given in the next section since the antidiagonal solutions 
appear to be particular cases of the mixed solutions discussed below. 
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5 Mixed solutions of the reflection equation 

We now look for invertible solutions with terms both in the diagonal and in the antidiagonal part 



m+n 



K{u) = {Hu)Eii + £i(u)E^ (5.1) 

i=l 

with at least one non zero 

Lemma 5.1 The case of so(2) being excluded, the solutions of the reflection equation of the form 
\5.1)) are all constant (up to a global normalisation function) . 

In the special case of so{2), the set of solutions of the reflection equation is the union of function- 
valued diagonal matrices and function-valued antidiagonal matrices. 

Proof: Projecting the reflection equations on the elementary matrices E pq ® E rs , one gets 

on Eij <g> Efl : (h{u) i^v) + k T (v) i^u)) + (kj(u) i t {v) - k 3 {v) i^u)) = (5.2) 

u + v V / u — v \ J 

on ® Eji : - tl£ (kj(u) £ s (v) + kj{v) £ 3 (u) \ - (h^u) £^ v ) _ k^v) £ 3 (uj) = 

(5.3) 

(-l)W / 

on E^ £g> Eji : — — — \ki{u) ki(v) + £i(u) £i{v) — kj{u) kj{v) — £j(u) £j(v] 



kAu) kAv) - kAu) kAv)) = (5.4) 

u — V v / 

on E^ ^Ej,: - ^(u) - ^(t») ^(u)) = (5.5) 

on £y <g> ^ : - 1 J OiOjOo \ (-l) [ %(u) £j(v) - (-l) b] £j(u) £&) 



««)«•)-««)«•)) =0 (5.6) 



u — V 

on ® ^ : (-lpl+M+bl^. / — fk^u) kj(v) - k T (vj) + 

+ (fei(u) fct(u) - fei(u) fcj(^) + 4(^) - (-l)^j(w) ^) 

- — JL - - (— (**(«) Hv) - kj(u) h(v)) + StlK ^ (h(v) - kj(v))) 

U + V + K \U — V v U — V + K ') 





h ( u ) h ( v ) _ k ( U ) fc (v ) + ^ _ £j(u) £j(v) =0 
(w + f )(it — v + k) x ' 1 



(5.7) 
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on Eij ® E- %i : (-lpl+M+blfl.fl. | ^ + ^ + ^ ^( u ) ^) - (_l)b1^ u ) ^(^)) + 

(kjiu) £j(v) + {-lf\ij{u) kj(v) + (kj(u) £j(v) + £j(u) kj(v)) 

(h{u) £j{v) - £j(u) h(v)) - StrK ^) £ _ {v) \ \ = (5 8) 



M + f + KV-U — t> M — V + 

U + f + K 

+ — (h(u) £i(v) + {-\f\iiiu) h(v) + (k- t (u) £i{v) + £ t (u) h{v)) 

U — V + K V U + V 



on £ y ® E i} : (-lpl+M+M^. ( 1 (k 3 (u) - + 



°_ kj (v) - kj(u) £ t (v)) + StrK(M) ^)) I = (5.9) 



U + V + K \U — V U — V + 



on 1®E U : -5—3 Ui(u) 40) - £- t {u) £i{v)) = (5.10) 

U — V V / 



on 1® S^: -= r (^(u) ^(u) + 4(u) - k T {u) l^v) - l^u) h{v) ) = (5.11 



W 2 — t> 2 



Of course, when several indices i,j coincide, the corresponding equations merge into a single one. 
Consider now a couple of indices such that i, J are all different. Then eq. (J5.5j) implies 

£i{u)£j{v) = £j(u)£ t (v), (5.12) 

the solution of which is given by £i(u) = £j(0) £(u) where £(u) is an arbitrary function, which can be 
factorised by a change of normalisation of the diagonal elements. Hence one can restrict K in (|5.1j) 
to have only constant antidiagonal elements without loss of generality. 

Suppose now that 4^0 for some i (purely diagonal solutions are not considered in this section). 
The general solution to the system formed by eqs. (J5.2j) and (J5.ll)) is then given by 2 

ki(u) = 7U + fcj(0) ; kj{u) = — fcj(0) and kj(u) = — ju + kj(0) (j ^ i,i) (5.13) 

where 7 is a constant. The inspection of the E w ® E^ coefficient leads to 7 = 0. 

It follows that all diagonal terms are constant. Hence the mixed solutions for the reflection matrix 

are necessarily constant. 

The case of so(2) is solved by direct calculation. ■ 
We now specify the exact form of these solutions. We establish (the case of so(2) being excluded): 

Proposition 5.2 Mixed solutions of the reflection equation for osp(m\n) exist only when m is even. 
They fall into two classes: 

CI: The so{m) block is diagonal. 



2 Actually, we have restricted ourselves to meromorphic functions on C. 
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The solutions are parametrised by n complex parameters. The matrix K is given by: 

Jvi — — X kj — 1 



forie{l,...,f}: 
for i G {m + 1, . . . , m + n} : 
C2: The sp(n) block is diagonal. 



4 = 

ki = sin(cti) 



with 



at + cti = 



(5.14) 



4 = e ft cos(ai) Pi + fit = 



T/ie solutions are parametrised by a couple of positive or null integers ni\ > m 2 such that 

m n n 

m\ + m 2 < — — 1 ; mi — m 2 < — and m\ — m 2 = — [mo<i 2] 

and by m — 2(mi + m 2 ) complex parameters. Setting ri\ = (n + 2m\ — 2m 2 )/4 ; £/ie matrix K 
is given by: 

ki — Kt — 1 



so(m) part: < 



sp(n) part: 



fori G {l,...,mx} : 
/or i G {mi + 1, . . . , mi + m 2 } : 
/or i G {mi + m 2 + 1, . . . , : 

/or i G {m + 1, . . . , m + rii} : 
for i G {m + n\ + 1, . . . , m + f } : 



4 = 4 = 

/u^ — /u^ — 1 

4 = 4 = 
— /l^ — 
4 = 4" 1 ,AeC 

nJfl — kj — X 

4 = 4 = 

k{ — k% — X 
4; = 4 = 



(5.15) 



(5.16) 



For so(m) (resp. sp(n) ), the mixed solutions are of the form C2 with n = and m 2 = mi (resp. CI 
with m = 0). 

Note that this classification is given up to a global normalisation function and up to a relabelling of 
the indices, as noticed in proposition ^. 1\ 

Proof: We have to consider only constant reflection matrices K. Then, extracting the residues at 
the poles in u, eqs. (15. 2j) to ()5.11|) reduce to (for i ^ i,j,J and j ^ j) 



4 + h) = o 


(5.17) 


k^ kj ~\~ ^L%^L% ^j^j ~~ ^ 


(5.18) 


hh - k 3 kj + #o((-l) [l] 44 - = 


(5.19) 


( (_i)W _ (-i)W)44 = o 


(5.20) 


/c^ A,^ k j k^ — 


(5.21) 


(^-(-1)^)4 = 


(5.22) 


(i + (-i)[% ) Mi = o 


(5.23) 


(i + (-i)[% ) fc,4 = o 


(5.24) 


(Jfei-fy) Strif = 


(5.25) 


4 Str AT = 


(5.26) 
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In the case of so(m) algebras, eq. (|5.22|) implies ki = kf for all 1 < i < m, from which (j5.21|) 
follows. Eq. ()5.17|) is a consequence of (|5.23|) and ()5.24|) . eq. ()5.25|) follows from ()5.26|) and (|5.19|) from 
(|5.18|) . Finally the mixed solutions in the so(m) case are characterised by the following equations: 

h = k- t ■ kdi = ; Tr K = (5.27) 
kf + = k) + £j£ 3 for all i, j (5.28) 

These constraints exclude the existence of mixed or antidiagonal solutions with odd m, since a matrix 
cannot be traceless if its diagonal has an odd number of non vanishing elements with equal squares. 

In the case of sp(n) algebras, eqs. (|5.23|) and (J5.24|) vanish, while eq. (|5.22j) implies fc« = —k^ for 
all 1 < i < n, from which (1537]) . (I5~2T1) . (f5~23J) and (l5~27)j) follow. Moreover, eqs. (|5~T%J) and flQ| 
are equal. Finally the mixed solutions in the sp(n) case are characterised by the following equations: 

ki = -ki ; TtK = (5.29) 
k\ + ^4 = k) + £j£j for all i, j (5.30) 

In the case of osp(l\n), one has to add the following equations (with 2 < i,t < 1 + n) 

k\ - kj - £ { £- % = (5.31) 

It follows from (j5T3T|l - (fP2jl and (PT^t that kf = for alH = 2, . . . , 1 + n. Since £ % Str K = 0, one 
has also k\ = 0. Therefore, there does not exist any mixed or antidiagonal (invertible) solution for 
osp(l\n). 

In the case of osp(m\n), (m > 1), eq. ()5.20|) shows that at least one antidiagonal part (orthogonal 
or symplectic) is zero: £i = for i £ [1, m] or l± — for i £ [m + 1, m + n]. This means that no pure 
antidiagonal solution exists for osp{m\n) superalgebras. We thus have to consider two cases. 
• Case 1 : the so(m) block is diagonal 

Eq. ()5.22|) with index j £ [m + l,m + n] implies ki = —ki for all i, which excludes the case of 
osp(m\n) with odd m. The remaining equations then lead to 

k] + £i£ T = k) + £j£ 3 for i,je[l,m + n] (5.33) 
kj = —kj for j £ [1, m + n) (5.34) 



• Case 2 : the sp(n) block is diagonal 

Eq. ()5.22|) with index j £ [l,m] implies ki = ki for all i. The remaining equations lead then to 

k\ + £i£ T = k) + £j£j for i, j £ [1, m + n] (5.35) 

fcj = fcj for j £ [1, m + n] (5.36) 

fe.^. = ki£- x = for i £ [1, m] (5.37) 

m m+n 

StriT=(^- )^ = (5.38) 

i=l i=m+l 
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Again, these constraints exclude the existence of mixed or antidiagonal solutions with odd m, since 
one cannot have a traceless matrix if its diagonal has an odd number of non vanishing elements with 
equal squares. 



For example, one finds the two following solutions for osp(4|2): 



/ 1 



1 



1 



k 5 4 

V 4 -k 5 ) 

where k\ + 44 = 1. For osp{2\4) the two solutions take the form 



\ 



and 



1 1 




4 




o 




1 






i 


V 


i 



/ 1 



-1 



\ 



h 4 
fc 4 4 

4 -fc 4 

4 -h J 



( 4 

£^ 



and 



-1 



1 / 



where kl + £ 3 £ e = 1 and k\ + 44 = 1- 



6 Analytical Bethe Ansatz for the so(n) and sp(n) open spin 
chains 

The main aim of this section is the derivation of the Bethe Ansatz equations for the so(n) and 
sp(n) iV-site open spin chains with non trivial reflection conditions by means of the analytical Bethe 
Ansatz method (see e.g. jSH E3 HHH OH E3 EH])- To construct the open chain transfer matrix we 
need to introduce the R matrix (371 EH1, which is a solution of the Yang-Baxter equation (|2.7|) . 
Normalisation of R matrix will be modified in order to connect it with the physical XXX type 
Hamiltonians considered hereafter; in addition, we set A = iu. We focus on the so(n) or sp(n) 
invariant R matrix given by: 

R(X) = A(A + ik)I + i(X + ik)P - i\Q. (6.1) 

The R matrix satisfies also crossing and unitarity, namely 

#i 2 (A)i? 12 (-A) = (A 2 + k 2 )(A 2 + 1)I, R 12 (X) =R[\(-X-tK). (6.2) 

where * is the transposition defined as in (|2.Hj) . 
The open chain transfer matrix is defined by [2] 

t(X) = Tr K+(X) T (A) K (X) f (A) , (6.3) 



II 



where K (A) is any solution of eq. and Kq(X) is a solution of a closely related reflection 

equation defined to be: 

R 12 (v - u) K{ x (u) R 12 (-u -v- 2in) K l 2 2 (v) = K l 2 2 (v) R 12 (-u -v- 2m) K£(u) Rn(v - u) (6.4) 

Tr denotes trace over the auxiliary space 0, and 

T (A) = R 0N (X)R 0N -i(X) ■ • ■ i? 02 (A)i? i(A) , T (A) = R 10 (X)R 20 (X) ■ • ■ R N ^ (X)R N0 (X) (6.5) 

It is clear that any solution K~(X, £ _ ) of where £~ are arbitrary boundary parameters, give 
rise to a solution K + (X) of ()6.4j) defined by K + (X) = K~(—X — m)', also depending on arbitrary 
boundary parameters £ + = 

To determine the eigenvalues of the transfer matrix and the corresponding Bethe Ansatz equa- 
tions, we employ the analytical Bethe Ansatz method [HTJIII2]- Namely, we impose certain constraints 
on the eigenvalues by exploiting the crossing symmetry of the model, the symmetry of the transfer 
matrix, the analyticity of the eigenvalues, and the fusion procedure for open spin chains. This then 
allows us to determine the eigenvalues by solving a set of coupled non-linear consistency equations 
or Bethe Ansatz equations. 

We first focus on the case with trivial boundaries, namely K~(X) = K + (X) = 1. We will in a second 
step derive the Bethe Ansatz equations for all diagonal solutions found in the previous sections, 
namely Dl, D2, D3, D4 and solve them in the thermodynamical limit. 

To obtain the necessary constraints, we recall that the fusion procedure for the open spin chain 
[SHI EH IHE] yields the fused transfer matrix 

t(X) = C(2A + 2%k) t(X) t(X + ik) - (C(A + iK)) 2N q(2X + m)g(-2A - 3m) , (6.6) 

where we define 

C(A) = (A + m)(A + i)(A-m)(A-i), (6.7) 

f(A-0(A-«) for so(n) 
I (A + X + in) for sp(n) 

In addition, from the crossing symmetry of the R matrix ()6.2|) it follows that (when K~ = K + = 1): 

t(X) =t(-X-iK). (6.9) 

The transfer matrix with K~ = K + = 1 is obviously so(n) (resp. sp(n)) invariant since the corre- 
sponding R matrix (|6.1j) is so(n) (resp. sp(n)) invariant. The symmetry of the transfer matrix makes 
the computation of its asymptotic behaviour (A — > oo) a relatively easy task. Finally, to implement 
the analyticity of the eigenvalues, we require that all poles must vanish. These constraints shall 
uniquely fix the eigenvalues. This is the basic outline of the analytical Bethe Ansatz method. 
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6.1 Eigenvalue of the pseudo- vacuum 



In order to compute the general eigenvalues we need to first define a reference state or "pseudo- 
vacuum". After finding its pseudo-energy eigenvalue, we will be able, with the help of the above 
discussed constraints, to derive the general eigenvalues of the transfer matrix. The pseudo-vacuum, 
which is an exact eigenstate of the transfer matrix, is the state with all "spins" up, i.e. 



N 



where |+) is the n-dimensional column vector 



1+) 



l+>< 



(6.10) 



i=l 



/ 1 \ 





(6.11) 



The action of the R matrix on the state |+) from the left gives rise to upper triangular matrices, 
whereas the action from the right on (+| = |+)' gives lower triangular matrices. It is obvious then 
that the corresponding action of the T (resp. T) on \u+) (resp. (o>+|) will also give rise to upper (resp. 
lower) triangular matrices. It is relatively easy, after some tedious algebra (see jHHl for a detailed 
example), to determine the action of the transfer matrix, t(X)\uj + ) = A°(A)|o; + ), where A°(A) is given 
by the following expression 



71-2 



A°(A) = a(X) 2N g (X)+b(X) 2N Y,9iW+cW 2N 9n-i(X) 



i=i 



with 



a(X) = (X + i)(X + iK), b(X) = X(X + iK), c(A) = A(A +in — i) = a(—X - in) 



(6.12) 



(6.13) 



m (A+f±jXA+jK) 
9oW = — /x , , n — ' 9n-l[X) 



(A + f t|)A 



g {-X - in) 



(6.14) 



(A + f)(A + f) ' w (A + f)(A + m-f) 

In g and g n -\ the upper sign corresponds to so(n) and the lower sign to spin). Moreover (when 
n = 2k), 



9iW 
9iW 



A(A+f ±|)(A + i/Q 
(A + f)(A + f)(A + 



< I < k, 



g n -i-i(—X — in), k < I < n — 1 (6.15) 

This expression is valid for so(2k) with the upper sign and for sp(2k) with the lower sign. For the 
so(2k + 1) chain the expressions for gi, I ^ k, are the same as in the so(2k) chain, except for 

A(A + in) 



9k{X) 



(A + f)(A + ^; 



9k{-X 



IK) 



(6.16) 
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6.2 Dressing functions 

Now that we have the expression for the pseudo-vacuum eigenvalue, in accordance with the general 
analytical Bethe Ansatz procedure, we make the following assumption for the structure of the general 
eigenvalues: 

n-2 

A(A) = a(X) 2N g (X)A (X) + b(X) 2N J>(A)A(A) + c(A) 2 V-i(A)A n _ 1 (A) (6.17) 

i=l 

where Ai(X), the so-called "dressing functions", will now be determined. 
We immediately get from the crossing symmetry of the transfer matrix (|6.9|) : 

Ai(X) = A n _ z _!(-A - m) 0<l<n-l. (6.18) 

Moreover, from the fusion relation (J6.6j) . we obtain the following identity by a comparison of the 
forms ()6.17|) for the initial and fused auxiliary spaces: 

A (A + zk)^-i(A) = 1. (6.19) 

Gathering the above two equations ()6.18|) . ()6.19|) we conclude 

Ao(A)A (-A) = 1 . (6.20) 

Additional constraints are obtained on the dressing functions from analyticity properties. Studying 
carefully the common poles of successive g^s, we deduce from the form of the gi functions ()6.14|) - 
f!6.16|) that gi and <#_i have common poles at A = — ^, therefore from analyticity requirements 

A l (~) = A i _ 1 {- t ±), i = l fc — 1. (6.21) 

The last relation is valid for so(n) and spin) as well. However, there is one extra constraint specific 
to so(2k + 1), namely 

A k {--) = A k . l {--). (6.22) 

Having deduced the necessary constraints for the dressing functions, we now determine them ex- 
plicitly. The dressing functions A\ are essentially characterised by a set of parameters {A^ | j = 
1, . . . , M^}, where the integer numbers are related to the eigenvalues of diagonal generators of 
so(n), sp(n). Defining these generators as: 

N 

S m = J2sf\ s^=Eu-E tt . (6.23) 

i=l 

The precise identification of M follows from the symmetry of the transfer matrix (see also [H5)> i n 
particular 

s (i) = M (i-i) _ M (0 ) / = !,..., k -2 (6.24) 



14 



which is valid for so{2k + 1), so(2k) and sp(2k). The two remaining quantum numbers are given by 



For so(2k + 1) : = M {k ~ 2) - M {k ~ l) ; S {k) = M {k ~ l) - M {k) ; 

For so{2k) : = M (fc ~ 2) - M (+) - M+ } ; S (k) = M (+) - M+ } ; 

For sp{2k) : = M (fc " 2) - M (fc - 1} ; S (k) = M (fc - 1} - 2M (fc) . (6.25) 

We have also = N (valid for both sp{n) and so{n) algebras). 

Let us point out that, away from the boundaries, the behavior of the chain is considered to be as 
in the bulk. Therefore, the above quantum numbers describe accurately enough the states of the 
system. The dressing functions are given by the following expressions: 

A. so(2k+ 1) 



^o(A) = I] 



3 2 " "3 2 



j=i A + Af ^ + § A - Af ) + \ 

A(A) = T] 1 m \,t * 

M A + Af + f A-Af + f 



M(i+1) a + A (m) + - - ^ A - A (m) iJ 



x 



n _j ^2 2 "j 1 2 2 7 -1 

j=l ^ + + g + 2 A ~~ A j +2 + 2 



^A + Af + f -i A-A^ + f -iA + Af } + f + i A - Af } + f + f , 
4t(A) = [ ^r-^ -^r^ m ~ m " 6-26 

11 X , \( k ) j_ik \ , ife \ , , ifc i \ \(*0 , ifc i v 7 

j = x at Aj i — 2" ^ — i — 2" i — 2 — 2 — i ~^ — 2 — 2 

together with A(A) = A n _/_i(— A — i/c) for / > /c. We recall that for so{2k + 1), k — k — \. 
B. so(2k) 

The dressing functions Ai(X) are the same as in the so{2k + 1) case for I — 0, . . . , k — 3, while 

A k -2{\) = n 

m(+) 

M(") 



A + Af - 2) + f 




A 


_ A j + T 


A + Af ~ 2) + f - 


- i 


A- 


j + 2" ~ * 


A + Aj +) + f - 


3i 
2 


A- 


\ (+) i ifc 3i 

A i + T ~ T 


A + Aj +) + f - 


i 
2 


A- 


\(+) i ife i 

A i + T ~~ 2 


A + Af } + f - 


3i 
2 


A- 


\ ( — ) i ife 3i 

A i + T ~ T 


A + Af } + f - 


t 
2 


A- 


A i + T - 2 
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A fc _i(A) 



A + Al +) + ik 



M(+) 

n A+A <+> 



2 _ T A ~~ \ 



(+) i ik 3i 

2 2 



2 ~ 2 A 



A j + T ~~ 2 



\ -\<f } + '4 + i A-A^' + g + j 



2 



7=1 A + A 



(-) 



T ~~ 2 A 



Ai- ' + ^ 



(6.27) 



with still ^4z(A) = A n _/_i(— A — z/c) for I > k — 1, with k, = k — 1. 



C. sp(2k) 

Similarly, the dressing functions are the same as in the so(2k + 1) case for I = 0, . . . , k — 3, and: 

M( fc - 2 ) , x(fc-2) , ifc \ \(*-2) ,ik 

A k . 2 (x) = n ' T y T 



; , A + Af- 2) + f -i A-Af- 2) + f -i 

x TT _ j + 2~ ~ ¥ A ~ A j + T ~ I 

j=1 A + Aj + Y - 5 A - Aj + y - 2 

4 m TT A + A i + T + 2 A ~ A j + T + 2 

^k-H/V J_J_ (fc-1) ifc j .(fc-l) ik 

j=l A + A j ~2 ~ 2 j + ~2 ~ 2 

x TT j +T~T A ~ A j + T~T /g 2gl 

3=1 A + A i + 2 + 2 A - A i + 2" + 2 

In addition Ai(X) = A n _i_i(— A — m) for £ > k — 1; note that now « = /c + 1. 



6.3 Bethe Ansatz equations for K =1 

The above dressing functions Ai satisfy all the imposed constraints and they are unambiguously 
defined. 

Requiring now the analyticity of the eigenvalues, we deduce the Bethe Ansatz equations. More specif- 
ically, successive A{s have common poles, which must disappear. Hence, the sum of corresponding 
residues of Ai and Ai + i in the eigenvalue expression ()6.17j) must be zero. The Bethe Ansatz equations 
immediately follow from this condition. 
Let us define 

A + — 

e,(A) = (6.29) 
A ~ ¥ 

for any x. Then, the Bethe Ansatz equations read: 
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A. so(2k + 1) 

«.(*!")"= II ^(Af'-Af'je^AW+Af) J] «-i(\ W - ) «-i(\ W + ) . 
i=ij¥* 3=1 

1= II e 2 (A?-A?)e 2 (A? + A?) J] II M*! " ^ + f) 

j=l,j^i T)=±l j=l 

l = 2,...,k-l 

i= n ^(Af-Af^e^Af+Af) n e -i(^ fc) - A r i) ) e -i( A ? ) + A r i) ) 

(6.30) 

B. so(2k) 

The first k — 3 equations are the same as in so(2k + 1), see eq. (|6.3Uj) . but the last three equations 
are modified, namely 

i = n ^t 2) - ^r 2) ) e 2 (\t 2) + xf ^ 2) ) n e -^ (fc ~ 2) - A f _3) ) e -^ (fe ~ 2) + a ? ~ 3) ) 

i=ij¥* 3=1 
Af(+) m(-) 
x n e -i(A; fe-2) - A<+>) C _ l( Af- 2 ' + A<+>) rj -i(AS'- 2 ' - A'"') e_ l( Af- 2 » + A<"») 

i= n e 2 (Ar-Af) e2 (Af)+Af) n M^-^vrf+ru^-i 6 - 31 ) 

3=ljV* 3=1 

C. sp(2k) 

The first /c — 2 equations are the same as in the so(2k + 1), while the last two equations are given by 
1 = II e ^ X t 1] - Xf- 1] ) e 2 (\t 1] + Xf -») J] e-^t 1] ~ Xf' 2) ) e^A?^ + Af' 2) ) 

3=ljV* 3=1 

x II e ^( A ? _1) - A ? } ) ^(Af-^+Af) 

3=1 

1= II ^(Af-Af^e^Af+Af) IlM^-^VrfH^)- (6-32) 

3 = ljV* 3=1 

For all these cases, we recover the rational limits of the equations given in jSUESj- 
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6.4 Eigenvalues and Bethe Ansatz equations for diagonal K 

Until now, we have considered the case of trivial boundary effects K + = K~ = 1. We here come 
to the main point of our derivation and insert non-trivial boundary effects. We shall then rederive 
modified Bethe Ansatz equations. We choose K~ to be one of the diagonal solutions Dl, D2, D3, D4. 
We consider, for simplicity but without loss of generality, K + = 1. Remark that the pseudo- vacuum 
remains an exact eigenstate after this modification. 

Let us rewrite the solutions Dl, D2, D3 and D4 in a slightly modified notation, which we are 
going to use from now on. 

Dl: The solution Dl can be written in the following form 

K(X) = diag(a, . . . , a, ft, . . . , ft) . (6.33) 

The number of a's is equal to the number of ft's, so that this solution exists only for the so(2k) and 
sp{2k) cases as stated in Proposition 13.11 and 

a(A) = -A + ^, /3(A) = A + i£, (6.34) 

where £ = - is the free boundary parameter. 
D2: Solution D2 can be written as 

K (A) = diag(a, ft, . . . , ft, 7) (6.35) 

with 

a(A) = ^±|i, ft(X) = 1, 7 (A) = ^±^, (6.36) 
A + if ! A + < n 

where fi = — f n = — y- are the boundary parameters which satisfy the constraint 



& + £n = «-l. (6.37) 

We remind that this solution exists for so(n) algebras, but not for sp(n). 
D3: Solution D3 has the form 

K(X) = diag(a, . . . , a, ft, . . . , (3, a, . . . , a). (6.38) 

The number of a's is 2m and the number of ft's is n — 2m, where the free integer parameter m obeys 
1 < m < k — 1 for so(2k), sp(2k) and 1 < m < k for so(2k + 1). Moreover 

a(A) = -A + if, /3(A)=A + if (6.39) 

where f = j — m has a fixed value, unlike the su(n) case where the corresponding solution has a free 
boundary parameter (see [4^ I41j). 
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D4: Finally solution D4, which holds only for the so(4) case, can be written in the following form 

K(X) = diag(a, (3,^,6) (6.40) 

where 

a(A) = (-A + z£_)(-A + z£ + ), /?(A) = (A + z£_)(-A + z£ + ), 

7 (a) = (-A+ie_)(A+ie+), *(A) = (A+ie-)(A+»e + ). (e.4i) 

= £ + = i are both free parameters. 

We now come to the explicit expression of the eigenvalues when K~ is one of the above mentioned 
solutions. We should point out that the dressing functions are related to the bulk behavior of the chain 
and thus they are form-invariant under changes of boundary conditions. Indeed what is modified in 
the expression of the eigenvalues (j6.17j) are the g% functions which characterise the boundary effects. 
We call the new g\ functions g\. 

Dl: As already mentioned, the solution Dl can only be applied in so(2k) and sp(2k). In this 
case we have 

<#(A) = (-A + z£)<?z(A), l = 0,...,k-l 

&(A) = (X + i^ + m)gi(X), I — k, ... ,2k — 1 (6.42) 

where g(X) are given by (|6.14p — ()6.16|) . The system with such boundaries has a residual symmetry 
sl(k) in both cases, which immediately follows from the structure of the corresponding K matrix. 

D2: We have 

~ m (-A + ifi) m ~ m (A + z^i + z) (A + z^i + z'k) 

(A + z£i) (A + z^i) (-A - ik + z£x + z) 

5/ (A) = (A alV^ ( A )' i = l,---.n-2 (6.43) 
(A + zCi) 

Again the gi(X) are given by ()6.14|) - ()6.16|) . Similarly, from the structure of the K matrix we conclude 
that the residual symmetry is so(n — 2) <g> so(2). 

D3: For the D3 solution we find the following modified g functions 

9i(X) = {-X + iOgi{X), Z = 0,...,m-1 

gx(X) = (A + — ± -)gi(X), l = m,...,n-m-l 

(A + - i ') 

9i(X) = (-X-iK-i^ jT— t k K fl(A), Z=n-m,...,n-l. (6.44) 

(A + y =F 2 J 

The symmetry of the transfer matrix for this K matrix is so(n — 2m) ® so(2m), (resp. sp(n — 2m) (g) 
sp(2m)). 
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D4: Finally for solution D4 the modified g functions are given by 

g (X) = (-A + i£_)(-A + i£ + )£ (A), ft (A) = (A + i£- + i)(-A + i£+)0i(A) 
£fe(A) = (A + i£+ + i)(-A + i£_)(7 2 (A), ft(A) = (A + ie-+i)(A + ze+ + ^ 3 (A). 

We now formulate the Bethe Ansatz equations for the general diagonal solutions. The only modifi- 
cations induced on equations ()6.30|) - (j6.32|) are the following for each solution: 

Dl The factor — e 2 ^ 1 +K (A) appears in the LHS of the k th Bethe equation. 

D2 The factor — e^ l+1 (A) appears in the LHS of the first Bethe equation. 

D3 The factor — e 2 ^ 1 +m (A) appears in the LHS of the m th Bethe equation with m = 1, . . . , k — 1 for 
so(2k + 1), sp(2k) and m = 1, . . . , k — 2 for so(2k). For so(2k), when m — k — 1, the factor 
— e ] ^ 1 (A) appears in the LHS of the (k — l) th and k th Bethe Ansatz equations. 

We treat solution D4 separately in the next section. 



7 Ground state and excitations 

The next step is to determine the ground state and the low-lying excitations of the model. One of 
the main aims of this work is indeed the computation of the scattering of the low-lying excitations off 
the boundaries. The bulk scattering for these models has been already studied in [22], nevertheless 
we are going to rederive these results as a check in the next section. 

We recall that the quantum numbers that describe a state are given by ()6.23|) . and the energy is 
derived via the relation H = ^-t(A) | a=o- It is given by 

e = -— y — -, — . (7.i) 



13 J 4 

In what follows we write the Bethe Ansatz equations for the ground state and the low-lying excitations 
(holes) of the models under study. Bethe Ansatz equations may in general only be solved in the 
thermodynamic limit N — > 00. In this limit, it is assumed that a state is described in particular by 
the density functions cr l (X) of the parameters xf . 

We here make the hypothesis that non trivial boundary effects do not modify the nature of the 
ground state and excited states but only the values of the Bethe parameters. 



A. so(2k+ 1) 

Let us first consider the so(2k + 1) case, for which the ground state consists of k — 1 filled Dirac seas, 
whereas the k th sea is filled with two-strings of the form Aq ± 4. Remark the shift |, instead of | 
usually expected, because the length of the k th root in the so(2k + 1) case is half the length of the 
other roots (see also [2Zj )• This leads us to rewrite the Bethe Ansatz equations for the ground state 



20 



and the low-lying excitations (holes) in the following form, inserting the two-strings contribution in 
the k th set: 

e l( Af>r« = - rT a 2 (A« - A<") e 2 (Af> + A?) J[ e_ l( A?> - Af ) e_,(A«" + Af) . 

e 1 (A!") = -n^ -Af)e 2 (AS" + A<") TT II «->(*?> - Af") e_ l( A?» + A^') 

j=l r=±l i=l 

(Z = 2,...,A:-2) 

Af( fe -i) 

ei(Af- 1} ) = - J] e 2 (AS fc - 1 )-A}*- 1) )e 2 (A i (fc - 1) + A5*- 1) ) 
i=i 

Af( fc - 2 ) 

3=1 

Af( fc ) 

x II M^""" - M^"" + *f Kl^"" - Af ) M (Af-" + Af ») 
i=i 

} ) = - II e ^ A ^ - A f ^ e ^ A ^ + A f ^ e ^ fc) " A f) ^ + A f } ) 
i=i 

x n «-i^ w - A r"» m< a . w + ^"V|(\ w - a?- i; ) e_ J( Af» + Af-") 

(7.2) 



For the general diagonal solutions of the reflection equation we have to multiply: 

- for D2, the LHS of the 1 st Bethe Ansatz equation with — e^ 1+1 ; 

- for D3, we multiply the LHS of the m th Bethe Ansatz equation with — e^l+rrn 1 < m < k — \. 

We are interested in the low-lying excitations, which are holes in the filled sea and are highest weight 
representations of so{2k + 1). We restrict ourselves here to the states with v"> holes in the I sea, 
which correspond to the vector representations of so{2k + 1), and also to holes in the k sea, which 
corresponds to the 2 fe -dimensional spinor representation (see also [HI])- We convert the sums into 
integrals by employing the following approximate relation jl2] 

1 M r°° 1 vW _ 1 1 

i=l Jo i=l 

where the correction terms take into account the holes located at values xf and the contribution 
at + . We shall denote by f(uj) the Fourier transform of any function /(A). 

Once we take the logarithm and the derivative of (|7.2)1 . we extract the densities from the equation 

£(u))a(u) = o(w) + -^(w) + ^G(u, (7.4) 
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id 

where a x (X) = — — lne x (A) and a x {uj) = e~~ . We have introduced 



a(A) 



( 2ai(A) \ 




(7(A) 



V / 

F(X), G(X,£) are k component vectors as well with 



/ * l (A) \ 
a l (X) 



(7.5) 



A') 



= a 1 (X)5 jl - tn(A) + a 2 (A) + ^ (a 2 (A - Xf ) + a 2 (A + Xf))<f y 

j'=i 

- £ ( ai (A - Xf ) + ai (A + Xf )) (<f ilI+1 + 5 jV _i) , (i = 1, . . . , A; - 2) 



F k -\X) = a 2 (A)-(ai(A) + a f (A))-^(a 1 (A-Af)+a 1 (A + Af))^ fc _ 2 

+ ( a 2( A - A ?) + a ^ x + - (( a ! + a i)( A - *?) + H + a i)( A + A f)) 5 ' 



/(0 



F fc (A) = 3ax(A) + a 2 (A) - (ai(A) + a f (A)) - ^ f(oi + a f )(A - Xf ) + (oi + a f )(A + 



+ ^ ((2ai + a 2 )(A - Xf ) + (2ai + a 2 )(A + Xf ))<$« 



(7.6) 



and the £ dependent part is 

for D2: 
for D3: 



G j (X,0 = -a^ 1+1 (X)S jlf 

G J (X,£,) = -a 2 ^ +m {X)5 jm (m=l,...,k-l). 



(7.7) 



Finally, the non vanishing entries of fC are 



= (1 + a 2 (u}))5ij - a 1 (u)(5 iij+1 + 5ij-i), i, j = 1, . . . , k - 1 
K.k-i,k(u) = £fc,fc-iM = -(fli(w) + a|(a>)), 
fckk{uj) = 1 + 2ai(w) + a 2 {uj). 



The solution of (|7.4|) has the form 



(7(A) = 2 e (A) + 1$ (A) + 1$ 1 (A,0 



(7.8) 



(7.9) 
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where e and $0,1 are k component vectors with 

k k 

?(u) = n il (u)a 1 (u>), <I*H = ^^H^H, 6»(w,0 = X)^v(a;)^'(w,e) (7.10) 

3=1 3=1 

1Z = JC^ 1 and e 7 is the energy of a hole in the j sea which can be written in terms of hyperbolic 
functions 



coah(fc \ j) ; 7.M 

cosh(A;-i)| ' J ' ' ' v ; 2cosh(A;-i)| 

sinhmin(i, cosh(/c — ~ — max(i,j))~ 



7e iJ -(w) = ea v 1J ' 2 - - — — i2^i±± ij = i,...,k-i 

JK ' cosh(A;- |)f sinhf ' ,J ' ' 



c 2 sinli ^— 

K 3h {u) = K kj {u) = T cosh(fe _ |) 2 fsinhr i = 1, 1 



sinh ^ 

Y 2 cosh f cosh(A;- |)| sinh f 



= — ., , , 2 u u ,.,, M (7-12) 



B. so(2k) 



In this case, the ground state consists of k filled Dirac seas of real strings. Therefore the Bethe 
Ansatz equations have exactly the same form as in ([6.31)1 . For the general diagonal solutions of the 
reflection equation, we have to multiply: 

- for Dl, the LHS of the k th Bethe Ansatz equation with —e^ +K , 

- for D2, the LHS of the I s * equation with — e^ i+1 , 

- for D3, the LHS of the m th equation with — for m — 1, . . . , k — 2, while the LHS of the 
[k — l) th and k th equations for m = k — 1 are multiplied by — e± . 

In this case as well, we restrict ourselves to states with i/W holes in the I sea. Note that now the 
spinor representation splits into two spinor representations of dimension 2 k ~ 1 (see also [3T]), and the 
holes in the +, — sea correspond exactly to these two spinor representations. The densities satisfy 
the same equation ()7.4|) as in the so(2k + 1) case with a and a given by ()7.5|) and 



coshfA; - 1 - 7')- 

g( "> = Jh ( *-i)g • ' 1 ' ^" ) = 2c M h(t-ir (7 ' 13) 



23 



i,0) 

F'(A) = a a (A)<y il -a 1 (A) + a 2 (A) + ^(a 2 (A-Af) + a 2 (A + Af))<y, j 

i=i 

- £ ( ai (A - Af ) + ai (A + Af )) (<f i)J+a + (i = 1, • • • , k - 3) 

F k - 2 (\) = a 2 (A)-2a!(A)+ J](a 2 (A-Af) + a 2 (A + Af))^_ 2 

i=i 

- ]T (a a (A - A<°) + a x (A + Xf )) (<V_ 3 + + <U 



i=i 

„(0 W (0 



F±(A) = 03(A) + £ (a 2 (A - Xf ) + a 2 (A + Xf ))<f, ± - £ ( fll (A - Xf ) + a x (A + Xf ))<^_ 2 

(7.14) 



j'=i i=i 



for Dl: G^A,£) = -a 2e+K (A)^ fe , 
for D2: G j (\,£) — -0( 2?1+1 )(A)^i, 

for D3- J G ^ A '^ = ~ a (2?+™)( A ) 5 i™ ( m= l,-..,fc-2), 

\ G J (A) = -a 1 (A)(^ 1 + 5 i , fc ) (m = fc-l). l ' ' 

Again, the non-vanishing entries of /C are 

JCij(u) = (1 + a 2 (u;))<% - ai(a;)(5ij + i + z, j = 1, . . . , fe - 2, 

/C__(cj) = JC ++ (u) = 1 + a 2 (w) and /C + _(cu) = = 0. (7.16) 

We solve equation ()7.4j) and find the densities a 1 which describe a Bethe Ansatz state. The solution 
of (J7.4)) has the same form as in (|7.9|) with 



sinhmin(i, j)^ cosh(A; — 1 -max(i,j))| 
cosh(/c — 1)^ sinh 2 



Kij{w) = e* — f-r: ^ . . .. 2 -, i,j = l,...,fc-2, 



, s ea sinh ^ 

^H=^M = T ooah(fc _ 1) | Biphr j = i,-.,*-2 

^ / n A / n sinh^f 



2 2 cosh | cosh(A; - 1) f sinh f 
el sinh(&-2)f 



K + » = = T 2cQshfcQsh(fc _ 1)!sinh| (7-17) 

Let us now consider the particular solution D4 for the so(4) case. The corresponding Bethe 
Ansatz equations, as in the bulk, are basically two copies of the XXX spin chain equations (so(4) = 
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su{2) ® su(2)), namely 

e^ + i(\ ( : ) )- 1 ei(\V) 2N+1 = U e ^ ) - X ^ e ^ ) + 4 ) ) ( 7 - 18 ) 

3=1 

where r = ±. It is obvious that the only representations that remain are the two two-dimensional 
spinor representations. The Bethe Ansatz equations are then two decoupled equations, as it is also 

evident from (j7.12|) : one has TZ + -(u) = 1Z- + {oj) = 0; moreover TZ ++ (u), 1Z (u;) are given by f)7. 17[) 

for k = 2, and 

F T {\) = ai(A) + a 2 (A) + ]T (a 2 (A - A] t) ) + a 2 (A + Af )) , G r (A, £ r ) = -a 2?T+1 (A). (7.19) 
C. sp(2k) 

The ground state in this case consists of k — 1 filled Dirac seas of two-strings (A^ ± |), and the k 
sea is filled with real strings. The Bethe Ansatz equations take the form 

e a (A? ] r +1 = ~ II ^(A! 1} - A«) e^(A« + A«) e 4 (A« - A« )e 4 (Af > + A«) 
i=i 

M( 2 ) 

x II -(Af - + *?) «-.(\ W - Af ) e_ 3 (Af» + Af >) , 

j=i 

e 2 (A?») = - II «5(\ W " Af) «K\ ffl + Af) e l( A» - A<")e 4 (A<<> + A?) 
»i=i 

x II II ^ - A " +T) ) ^(A? + A^) e _ 3 (A?) - Af^>) e- 3 (Af) + A^>) 

T=±l J=l 

l = 2,...,k-l 

MM 

e 2 (Af } ) =-II^* ) -^ ) ) ^(A^ + Af) 

3=1 

x II e-i^*' - X f e -i( A ?° + ^VsCA? - Af" 1 )) e_ 3 (Af } + A^) 
j'=i 

(7.20) 

where for the general diagonal solutions of the reflection equation we have to multiply: 

- for Dl, the LHS of the k th equation with — e"Tt +re , 

- for D3, the LHS of the m th equations with e 2 ^ 1 +m+1 e 2 ^ 1 +m _ 1 , m = 1, . . . , k — 1. 
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Again the densities for the state with holes in the / sea satisfy the same equations (|7.4|) as in the 
so{2k + 1) case with a given by (|7.5|) . a J (A) = 2a 2 (A)5,i and 



cosh(A; + 1 — j)\ 
2 cosh | cosh(fc + 1)| ; 



= o„„„uo,„„„u/>. I iw » J = 1,...,* (7.21, 



F J '(A) = 3a 2 (A)+a 4 (A)-2(a 1 (A)+a 3 (A)) + (a 1 (A) + a 3 (A))5 il 
+ ^ ((2a 2 + a 4 )(A - Af ) + (2a 2 + a 4 )(A + \f])8ji 



3=1 

z/(0 



- ^ ((ai + a 3 )(A - Af ) + ( ai + a 3 )(A + If )) + ^_x), 

3=1 

(j = 1) 

F fc (A) = a 2 (A) + a 4 (A)-(a 1 (A) + a 3 (A)) + ^(a 4 (A-Af) + a 4 (A + Af))5, 



/A: 



,(0 



- ^ (( fll + a 3 )(A - If) + ( fll + a 3 )(A + xf))^, (7.22) 

3=1 

for Dl: G3'(A,0 = -a2€+«(A)<y ifc , 

forD3: ©'(A, £) = -(a 2 ^ +m+ i(A) + a 2 £ +m -i(A))<5j m (m = 1, fc - 1) (7.23) 

The non-vanishing entries of /C are 

Kij{u) = (1 + a 2 (u;)) 2 % - (ax(a;) + a 3 (u;))(5 i j + i + 4j-i), i, j = 1, • • • , k - 1 
iC k -i,k(u) = K. k}k -i(uj) = -(ai(a?) + a 3 (a>)), 

/C fcfc H = l + a 4 H- (7.24) 

As in the previous cases the solution of (|7.4|) has the form (|7.9)L ()7. 10|) with 

- . , sinhminfz, ?')^ coshffc + 1 — max(i, , 

8 Scattering 

Having obtained the excitations of the model, we are ready to compute the complete boundary S 
matrix. For this purpose we follow the formulation developed by Korepin, and later by Andrei and 
Destri HI] • First we have to implement the so-called quantisation condition, 

{e 2iNpl S-l)\\f) = Q (8.1) 



26 



where p l is the momentum of the particle (in our case, the hole) with rapidity A^. For the case of v 
(even) holes in I sea we insert the integrated density ()7.9|) into the quantisation condition (|8.1|) . We 
use the dispersion relation 

= ~ p, < A > < 8 - 2 > 

and the sum rule iV / c?Acr(A) G Z+. The boundary 5 matrix for the right boundary is taken to be 
Jo 

one of the diagonal solutions Dl, D2, D3, D4, whereas the boundary S matrix for the left boundary is 
proportional to unit. We end up with the following expression for the boundary scattering amplitudes: 



a +l a~ l 



expfeirN J * d\(a l (\) -2e'(A))} (8.3) 
with 

a- l (\,0 = h(\)h(\,0, a +l (\)=k (X) (8.4) 

where or 1 is the first element of the diagonal boundary S matrix. It is obvious that a +l has no £ 
dependence and realises just the overall factor in front of the unit matrix at the left boundary (recall 
that K + = 1). Moreover, 

fco(Ai) =exp\iTrj ' d\& (\)}, fc^Ai, = exp {2m J ^AfKAU)} (8.5) 

with $[, 1 given by ([7. 10)1 . We finally restrict ourselves to I — 1 in the first sea and we write the latter 
expression in term of the Fourier transform (|7.1()j) of $g 1j 

fc (A) = exp{-i |°° ^(^e-^}, fc 1 (A,0 = exp{-^°°^<ll^,Oe^ A }. (8.6) 
In what follows we express the scattering amplitudes in terms of T-functions. 

A. so(n) 

Before we write down the explicit expressions for the boundary S matrices, let us recall the form of 
the exact bulk S matrix. It is easy to compute the scattering amplitude between two holes (vectors) 
in the first sea. The bulk scattering amplitude comes from the contribution of the terms of $J given 
by eqs. (|7.6|) . ([7.10)1 . ()7.12|) . (|7.14|) . with argument A±Aj. After some algebra and using the following 
identity 

i r du r(H±i) 

ln T^TsT ( 8 - 7 ) 



2 J to coshf r(^±5, 
we conclude that the hole-hole scattering amplitude is given by the expression 

tanvr(^) T(^) ^ + I) r(^±l) r(£§ + §) 



' S " (Aj tanTr(^) T(^) T(A, + k) T(^) r(=^ + ±) 

\ n— 2 / \n— 2/ vn— 2 2/ \ n— 2 / \ n— 2 2' 
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The explicit bulk S matrix then has the following structure 3 

s ( x ) = /-x S °w\ — tt(^(^ + «)I + ( iX + K ) p ~ iX Q) ( 8 - 9 ) 

(l\ + K)(l\ + 1) 

Now, we give the expressions for the boundary S matrix, which follow from (|8.5)L ()8.6|) . and the 
duplication formula for the T function 

2 2 "- 1 r(x + \)T{x) = 7r^r(2x). (8.10) 

The ^-independent part of the overall factor k , eq. (|8.5j) . is given by 

p/ a \ p/ -a _|_ 3x p/iX+l/2 , 3x p/ -iA+l/2 . 1\ 

/ _ V f\\ ^n-2> \n-2*&> \ n-2 '4/ L \ n-2 ' 2 > /o i i \ 

fc (A) - W)-— j- r( -iX+l/2 3 x r /iA+l/2 7 1\ {S - U) 

1 ^n-2^ 1 y n -2 ^ i> 1 \—^2 h V 1 I n-2 + 2~J 

where 

/iA+1/2 1\ • /iA— 1/2 , 1\ • / iA i 1\ 

Sin7T £ 7 sm7r 4~ + ^ sm7r + 7 

y- / \ \ v n— 2 4_^ V n-2 2> 'n- 2 A' / o -i c\\ 

1 0\ A ) — . /iA-1/2 , 1\ • /iA+1/2 l \ .i,,./ iA 1\ ■ \°' lZ ) 

Note that our solution includes the necessary CDD factors both for the bulk and boundary matrices 
(see also |2TJ|27|). The expression for the ^-dependant part ki depends on which solutions Dl, D2, 
D3, D4 we consider: 

Dl: 

fci(A, o = 1 n r* 2 1 ;: 2 , j (8.13) 

1 v n-2 ~ 2> L \ n-2 ~ / 

where £' = £ — | is the renormalised boundary parameter (see also |H]). We also compute the 
(5 element of the K matrix ()6.33|) by employing the "duality" transformation £ — > — £, and the 
symmetry of the matrix and of the transfer matrix (see also [2]). In particular, we obtain a set 
of Bethe Ansatz equations for £ — ► — £, which allows to determine the difference 

$1(^,-0- $1(^,0 = e- (2? - 1) ^ (8.14) 

and consequently 

r(A,o 



a-(A,0 



e^(A). (8.15) 



This provides us with a consistency check for our procedure since one obtains independently the 
exact ratio between different elements of the reflection matrix. Thus, we have completely determined 
the K matrix that corresponds to the solution Dl. 

3 We do not compute all the eigenvalues of the bulk S matrix via the Bethe Ansatz. Such a general computation, 
in the bulk, is beyond the scope of this work. 
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D2: 



*i(A,0 



tan7r 



n-2 



r(^# + ±) r(=^f + i) r(^±# + i) r( 



n-2 



n-2 



-»A+g; 

n-2 



tan7r 



■ n-2 ■ 



-»A+& 1. 
n-2 2' 



n-2 



jA+g; i 

n-2 ~ l ~ 2 



) r(^) 



(8.16) 



where £i = £i — | and ^ = £ ra + | are the renormalised boundary parameters. Remark that the 
constraint £i+£ n = k—1 is also true for the renormalised boundary parameters, namely £i+£„ = k—1. 
Similarly we employ the duality transformation which, for this solution, reads £i 
determine the difference 



,-(2fi-l)' 



+ e 



(2^+1) t 



£ n . We then 
(8.17) 



and the ratio 



7 (A,£l>£n) 

a~(A,fi,fn) 



e 2 ^(A)e_2«;(A). 



We have not been able to 



Again we have here an independent consistency check. 
Of course, we need to determine one further ratio, namely c , , , , ( . 
explicitly extract this information from the Bethe Ansatz formulation. Nevertheless, since the K 
matrix is a solution of the reflection equation, the only choice we have for this ratio is 



r(A,£i,e«) 



e 2^(A). 



8.19) 



a (A,fi,f n 

With this, we have completed the derivation of the K matrix that corresponds to the solution D2 
D3: The ^-dependent part of the K matrix overall factor that corresponds to D3 is 



fci(A,0 



r ( i\+e | 1\ p/ -iA+g' I i\ p/ -q+l/2 , 3\ p/ U+1/2 , 1\ 
1 ^ n-2 ^ 2> 1 \ n-2 ^ L ) 1 \ n-2 ^ D L \ n-2 4/ 



n-2 

1 \ n-2 2/ V n-2 



1^ 



n-2 
p/ iA+1/2 | 
1 v n-2 ~ l ~ 



3 - 
-) 



r( 



-iA+1/2 , lx 
n-2 4' 



(8.20) 



The total overall factor for this solution is 



fco(AMA,0 



>1(A) 



-iA ' 



:; , p pA+1/2 _|_ l- 



_L_ 

Vn-2 ~ 4 



r(^) r(^ + f) r( 



r(^+i) r 



X- 



n-2 



iA i 3\ 
4^ 

-iA+g' 
n-2 



-iX+1/2 , lx 
n-2 4' 



p/ -tA+l/2 
V n-2 



r 



-tA+l/2 
n-2 



+ 



(8.21) 



where 



n(A) 



/iA+l 

sm7r 

\ n— 



tA+l/2 
n-2 



T sm7r 



• q-i/2 

■ n-2 



+ |) sin^ + i 



/iA-1/2 . 1\ 



sin 7T 



•iA+1/2 

. n-2 



Sill 7T 



iA 
n-2 



(8.221 



The latter expression agrees with the one found in jzl], although we find slightly different CDD 
factors. Again £' = £ — ~ is the renormalised boundary parameter which now must have the fixed 
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value £' = | — m, to ensure the resulting .K" matrix satisfies the reflection equation. Unfortunately, 
in this case, we cannot employ a duality transformation to derive the ratio (|Q8j) . However, since 
the K matrix is a solution of the reflection equation the ratio must be 



r(A) 

a" (A) 



e 2 £'(A). 



(8.23) 



D4: For this solution the ^-independent part is given by (j8.11|) with n = 4, while the ^-dependent 
part is given by (r = ±): 



MA,e: 



2 'A' 

rf^ + |) 



(8.24) 



2 1 A' V 2 

We exploit the duality transformation one more time, £ T — > — £ r , to calculate the ratio 

/57(A) 



MA) 



e 2c;(A) 



(8.25) 



with the renormalised boundary parameter = ^ T — \. We recall that the two spinor representations 
are two-dimensional, and the corresponding K matrices are of course two-dimensional with 



K~(X,^) = diag(a;(X,C),p-(X,e T )) • 



(8.26) 



In other words, we have obtained two copies of the XXX boundary S matrix, with two free boundary 
parameters 

B. sp(n) 

The corresponding bulk scattering amplitude for the sp(n) is given by the following expression 



Sb(A) 



tan7r 



■ n+2 , 



tan7r( 



iA+l ' 
n+2 ■ 



T(— + -) Ti 

\n+2> Vn+2 ~ 2 > 



-iA+l ' 
n+2 - 



iA+l ' 
n+2 ■ 



p( iA+l , 1 
V n+2 2 
p/ -tA+1 _|_ 
V n+2 "'" 2 



IT ^(A) 



where 



5 6 (A) = exp 



00 rfw cosh ^ e~ iujX 



u 2 cosh f cosh ^ 2 



(8.27) 



(8.28) 



Again the structure of the S matrix is given by (J8.9)) with So given above. fco(A) is again given by the 
integral representation ()8.6|) using (|7.10|) . (|7.22|) and (|7.25|) . The corresponding ^-dependent parts 
for each solution are given respectively by: 



Dl: 



h(x,e) 



L V n+2 ~ 7.) L V 



-»A+g' 
n+2 



-iA+g' 
n+2 



1 V n+2 



1) 



(8.29) 
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where £' = £ — 1 is the renormalised boundary parameter. We also compute the f3 element of the K 
matrix by employing the duality transformation (£ — > — £) and we obtain the difference 

$K^-0- = e-«-=0* (8.30) 

and consequently 



e 2e (A). (8.3i; 



D3: 



tanTr(^-i) r(*g + §) r(^ + 1) r(^g£ + J) r ( ^ + i; 



i /\ pl\ _ — v n+2 A' v n+2 ' 2^ v n+2 ' / v n+2 ' 4^ v n+2 ' 4' /o oo\ 

11 " tan^ + i)r(^ + i) r(fijf + i) r ( ^ + i) + f) ' 1 ' 1 

= £ — 1 is the renormalised boundary parameter which must have the fixed value £' = | — m, 
in order for the resulting K matrix to satisfy the reflection equation. This result also agrees with 
the ^-dependent part of the corresponding result of [21] ■ Unfortunately we cannot employ a duality 
symmetry in this case to derive the ratio - (|6.38p . but again, since the K matrix is a solution of the 
reflection equation, the ratio must be 

^ = e^(A). (8.33) 
a (A) 



9 Conclusion 

We have established here the classification of diagonal, antidiagonal and mixed solutions to the 
reflection equations based on the Yangian R matrices for Lie (super) algebras so(m), spin) and 
osp(m\n). The next step would then be to obtain the complete classification of K matrices, still 
associated with one-dimensional quantum boundary space, which is technically more complicated. 
The cases of higher dimensional bulk representations, as well as operator valued i^-matrices remain 
to be classified. 

Some of the solutions we obtained have then been used to compute the spectrum and scattering 
datas for integrable spin chain systems with non trivial boundary conditions. Only the Lie algebra 
case has been considered here; the case of super Lie algebras osp(m\n) can now be envisioned: indeed, 
a suitable redefinition of indices allows us to consider a similar exact ferromagnetic pseudo-vacuum 
(see jH]), providing us with the starting point for the analytical Bethe Ansatz. 

Remark also that the analytical Bethe Ansatz construction proposed here remains valid for upper 
triangular reflection matrices. Such matrices can for instance be obtained by conjugating our diagonal 
solutions by triangular matrices (following lemma E3$ ■ 

In addition, the use of the analytical Bethe Ansatz method precluded applications to non-diagonal 
reflection matrices, for which no exact pseudo-vacuum state is available. However, this situation can 
be approached using the methods developed e.g. in jlHlllIj. For an alternative approach, see [48J. 
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